Nonlinear Integral Equations for 

Thermodynamics of the U q (sl(r + 1 
Perk-Schultz Model 



Zengo Tsuboi *and Minoru Takahashi ' 
Institute for Solid State Physics, University of Tokyo, 
Kashiwanoha 5-1-5, Kashiwa, Chiba, 277-8581, Japan 



Abstract 

We propose a system of nonlinear integral equations (NLIE) which 
describes the thermodynamics of the U q {sl{r + 1)) Perk-Schultz model. 
These NLIE correspond to a trigonometric analogue of our previous 
result, and contain only r unknown functions. In particular, they re- 
duce to Takahashi's NLIE for the XX Z spin chain if r = 1. We also 
calculate the high temperature expansion of the free energy. In par- 
ticular for r = 1 case, we have succeeded to derive the coefficients of 
order 0((^) 99 ). 
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1 Introduction 

Traditionally, thermodynamic Bethe ansatz (TBA) equations have been used 
to analyze thermodynamics of various kind of solvable lattice models pQ. 
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However TBA equations contain, in general, an infinite number of unknown 
functions, and are difficult to deal with. To overcome this difficulty, alterna- 
tive nonlinear integral equations (NLIE), which contain only a finite number 
of unknown functions, have been considered. For example, for XX Z-model 
(or XY Z-model), Klumper PUIS], Destri and de Vega [3] derived a NLIE, 
which contains two unknown functions. 

Recently, the second author derived another type of NLIE for the 
XX Z-model. The number of the unknown functions is only one, which 
corresponds to the rank of the underlying algebra U q (sl(2)). Moreover the 
same NLIE was rederived [S] from the T-system (a system of functional 
relations among transfer matrices) IE] of the quantum transfer matrix 
(QTM) from a fugacity expansion formula (THj. We find jTT] that this 
NLIE j^j is very useful to calculate the high temperature expansion of the 
free energy. In view of this situation, the first author derived NLIE for 
the osp(l\2s) model [12], the sl(r + 1) Uimin-Sutherland model [T3] and 
the spin | Heisenberg model [14] . The numbers of the unknown functions 
of these NLIE are s,r, 1, respectively, which correspond to the rank of the 
underlying algebras. In particular, the NLIE in refs. ^3] and ^3] contain 
the rational limit of Takahashi's NLIE [S] as special cases r = 1 or s = 1 
respectively. The high temperature expansion of the free energy from our 
NLIE [13] was applied [IH]-[II1] to spin ladder models, and good agreements 
were seen between experimental data and theoretical results. 

The Perk-Schultz model [201 [21] is a graded vertex model associated with 

the superalgebra U q (sl(r + l\s + 1)). In this paper, we consider a special 

case (s = —1) of it: the U q (sl(r + 1)) Perk-Schultz model [22]- This model 
is a natural multicomponent generalization of the XX Z model which is very 
important in physics. In fact the hamiltonian of the Perk-Schultz model 
(eqs. (]2.6[) and ([2.7]) ) reduces to the one for the XX Z model eq. ()4.2|) if 
r = 1. Moreover the rational limit of the Perk-Schultz model is the Uimin- 
Sutherland model [23123], which we considered previously [T3| . 

The purpose of this paper is to derive NLIE for U q (sl(r + 1)) Perk-Schultz 
model [271] l2*T] . This NLIE contains only r unknown functions, and corre- 
sponds to a trigonometric analogue of the NLIE in ref. ^H]- Moreover it 
reduces to Takahashi's NLIE [S] for XX Z-model when r = 1. 

The shape of TBA equations varies with q [23 EH] when q is a root of 
unity since the corresponding T-system (or F-system), from which the cor- 
responding TBA equations are derived, truncates with respect to a variable 
m (cf. eq. (12.20]) ). In contrast, our new NLIE eq. ([3.6]) will be also valid as 
it stand even when q is a root of unity since we use only the first part of the 
T-system (eq. (12.20)) for m — 1). 
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In 2, we introduce the U q (sl(r + 1)) Perk-Schultz model [213121] and define 
the QTM and the T-system. In 3, we derive our new NLIE eq. (|3.fij) taking 
note on the periodicity of the QTM. We will write only the outline of the 
derivation since it is similar to the ones in refs. [Bj and in 4, we will 
calculate the high temperature expansion of the free energy from our new 
NLIE. In particular, r = 1 case is a detailed explanation of a calculation in 
the letter [H], which has been extended to the XX Z-model case now. It will 
be difficult to derive the same results from the traditional TBA equations. 
Section 5 is devoted to concluding remarks. 



2 T-system and QTM Method 

We will introduce the quantum transfer matrix (QTM) method j^j, |2*7j - 

|22|, 0121 and the T-system for U q (sl{r + 1)) Perk-Schultz model. The 
QTM analyses of the Perk-Schultz model was done in ref. [3U] (see also, ref. 

ISUE2]). 

The i?-matrix of the model is given as 

r+1 r+1 r+1 r+1 

r w = EEEE R Z'» Eai ' a2 ® EblM i (2-i) 

»1=1 <J2 = 1 6i=l &2 = 1 

where E a,b isar + lbyr + 1 matrix whose element is given as (E a,b ) it j = 
SaiSbj] RZlfoiv) is defined as 

KIM = [v + 1]„ (2.2) 
K>) = [v} q (a^b), (2.3) 

<» = <l sisnia - b)v {a*b), (2.4) 

where v e C; a, b <E {1, 2, . . . , r + 1}; [v] q = (<f - q~ v )/(q - q' 1 ); q = e 71 . Let 
L be a positive integer (the number of lattice sites). The row-to- row transfer 
matrix on Cf +1 is defined as |BS] 

t(v) = tr (R 0L (v) ■ ■■R 02 (v)R 01 (v)). (2.5) 

The hamiltonian H = Hq + H c h of the model has two parts. The first part 
Hq is given by 

j=l a=l 
r+1 r+1 s 

E (sigKa-^sinhTyif^ (2.6) 

o=l 6 = 1 ^ 

a 7^ b 
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where we adopt periodic boundary condition E^ +1 = E^ b . The second part 
is the chemical potential term, 

L r+l 

H ch = -J2J2^ E T a - ( 2 - 7 ) 

j=l a=l 

Let R(v) be '90 degree rotation'of R(v): 

r+l r+l r+l r+l 

ai=l a2=l 6i=l 62=1 
r+l r+l r+l r+l 

ai=l 02=1 6i=l 62=1 
r+l r+l r+l r+l 

= EEEE ^;»£ 62 ' 61 ® ^ ai,a2 - (2-8) 

ai=l a2=l 6i=l 62=1 

Namely we have -R^fe^) = ^ll'al( v )- Let ^ be a positive even integer (the 
Trotter number). We define the QTM on Cf +1 as 

_H J _ _ 

tQTuyv) = trje t R NJ (u N + iv)R N ^ ltj (u N - iv) • • • 

R^j(u N + iv)R 3 j(u N — iv)R 2 ,j(u N + iv)Rij(u N — iv),(2.9) 

where j G {1,2, ...,L}; = — (T: temperature). The matrix ele- 

ments of the QTM are given as 

{a fc ,/3fe} 

N 

w<r.;S = E e ^ n k + w)^:^ («* - <«)(2.n) 

{v fc } fc=i 

where vn+i = v\ and Vk,ak,j3k £ {1, 2, . . . , r + 1}. We can express |H] the 
free energy per site in terms of only the largest eigenvalue Ai of the QTM 
eq. (EH) at v = 0: 

f = -T lim logA 1; (2.12) 

where the Boltzmann constant is set to 1. 

The eigenvalue formula of the row-to-row transfer matrix eq. (|2.5jl was 
derived |H4*| I2T] by the algebraic Bethe ansatz. On the other hand, the 
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eigenvalue formula T[ (v) of the QTM eq. ()2.9j) is conjectured to be given by 
replacing the vacuum part (the eigenvalue on the pseudo- vacuum |0 >) of the 
eigenvalue formula of the row-to-row transfer matrix with that of the QTM 
jHSJ- We adapt the following pseudo- vacuum: |0 >:= v r+1 <g> vi <g) • • • ® v r+1 <g> 
v 1 G Cf£, where v 1 :=' (1, 0, . . . , 0) G C r+1 and v r+1 :=' (0, . . . , 0, 1) G C r+1 . 
Indeed |0 > is an eigenvector for ^qtm^): 

r+1 

t QTM {v)\0>=J2Mv)\0>, (2.13) 

a=l 

where the eigenvalue consists of the functions: 

ipa(v) = e^(f) + (v + id a)r+1 )<f)-(v -id ajl ) for a G {1, 2, . . . , r + 1},(2.14) 

'smr](v ± zun) 



sinh?7 



Then the general eigenvalue formula T^ 1 "* (t>) of the QTM eq. ()2.9|) will be 
(cf. refs. EH- El) 
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Ti 1 \v) = J2^d;v). (2.15) 

d=l 

Here the functions {z(d;v)} are defined as 

i 

z(a;v) = tp a (v)- 



t Q a -i(v-i(a+l))Q a (v-i(a-2)) 



Qa-i{v - |(a - l))Q a (v - fa) 

for a G {1,2,..., r + 1}, (2.16) 

where Q a (v) = U^i^v(v - v^); M a G Z> ; Q (v) = Q r+ i(v) = 1. 
is a solution of the Bethe ansatz equation (BAE) 

Mv { : ] + - 2 a) _ Q a -i(v i : ) +i)Q a (v { : ) -i)Q a+ i(v^+^ 



^ a+1 (4 a) + §a) Q a ^(v { k a) - \ )g a (4 a) + *)Q a+1 (4 a) - j) 
for k <E {1, 2, ... , M a } and a G {1, 2, . . . , r}. 

Some remarks on eq. ()2.15|) are in order. For r = 1 case, eq. (|2.15j) can be 
proved jSHl by algebraic Bethe ansatz. As for r > 1 case, we have checked 
that eq. ()2.15|) agrees with the numerical diagonalization of eq. (J2.9J) for 
r = 2 and small Trotter numbers N. In ref. j3D] the algebraic Bethe ansatz 
was executed for one particle state, from which an eigenvalue formula of 
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the QTM for the Perk-Schultz model was conjectured, although a different 
pseudo-vacuum was adopted. 

As an auxiliary function, we define (cf. Bazhanov-Reshetikhin formula in 
ref. EH!): 



T ln\v) = \{\{z{d hk ;v--{a-m-2j + 2k)), (2.18) 

where the summation is taken over dj t k G {l,2,...,r + l} such that dj )k ~< 
dj + x tk and d jtk < d jtk +i (1 -< 2 -< • • • -< r + 1); m G Z> x ; a G {1,2, .. . ,r}. 
This function contains T^\v) eq. ()2.9|) as a special case (a, to) = (1, 1). It 
is related to a fusion [3H] hierarchy of the QTM. We can show that the poles 
of Tm\v) are spurious under the BAE eq. (|2.17|) . For a G {1, 2, . . . , r} and 
m G Z>i, we shall normalize eq. ()2.18|) as Tm\v) = Tm\v) j ' Mm\v) , where 



x n n * > ^ — o^+c* y 1 — (2.i9) 

3=1 k=l 

One can show that Tm\v) satisfies so called T-system [E] 

T^(v + l -)f£\v -'-)= T^Uv^Uv) + f^\v)f^\v), (2.20) 



for a G {1,2, ... ,r} and m G Z> 



where 



T (a) (w) = l for aeZ>i, 

~ (0) 0_(^ + f^ + (^-fO 

T^ +l \v) = e t for m G Z> x . 



3 Nonlinear Integral Equations with Only a 
Finite Number of Unknown Functions 

It is known jH] that Takahashi's NLIE [S], which describes thermodynamics of 
the AAZ-model, can be derived from the T-system of the QTM. The number 
of unknown functions for this NLIE is only one, which corresponds to the 
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(a,h)=(1,0) 
a,h = 2,0 
a,h = 1,3 
a,h = 2,3 
a,h = 1,6 
(a,h)=(2,6) 



Figure 1: Location of the roots {v^} 



1,2) of the BAE for U q (sl(3)) 



(r] = 2.5, N = 12, un = —0.05, /ii = h, jj, 2 = 0, /i3 = — h), which will give the 
largest eigenvalue of the QTM ^qtm^) at t> = 0. Color 1 roots {?4 } and 
color 2 roots {v k } form six one-strings. Only the roots in the fundamental 
domain Ret> G [— j-, f-] are exhibited. See also Fig. 1 in ref. [T3| for r\ = 0, 
/i = case. 



rank of the underlying algebra U q (sl(2j). In this section, we will derive the 

NLIE eq. f!3.6p for U q (sl(r + 1)), which contain only r unknown functions, 
from the T-system eq. ()2.20|) . 

From a numerical analysis for finite N,ujy,r, we conjecture that a one- 
string solution (for every color) in the sector y = Mi = M 2 = ■ ■ ■ = M r 
of the BAE eq. ()2.17|) gives the largest eigenvalue of the QTM eq. ()2.9|) at 
v = 0. From now on, we will consider only this one-string solution. Thus 
7f } (0) should give the largest eigenvalue A x . We expect that the following 
conjecture is valid for this one-string solution (cf. Figs. dandEJ). 



Conjecture 3.1. For small (\u^\ <C I), all the zeros {z[ } ofT[ c '[v) 
(a G {1, 2, . . . , r}) are located outside of a domain (the physical strip) Imv G 
[— n, o]- To be precise, they are located near the lines Imv 



n(a), 



_|_l+a 
- 1 - 2 • 
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(a,h) = (1,0) 
a,h = (2,0) 
a,h = 1,3 
a,h = (2,3) 
a,h = 1,6 
(a,h = (2,6) 



Figure 2: Location of the zeros of T{ 



v a 



1,2) corresponding to 



roots in Fig. Only the zeros in the fundamental domain Rev G [—j-, 
are exhibited. See also Figs. 2 and 3 in ref. ^3] for rj — 0, h — case. 



S 



We expect that this conjecture is also valid in the Trotter limit iV — > oo. 
f£\v) has poles only at ±j3$ n : (3$ n = ^i + iu N + ^ (n G Z) whose order 
is at most N/2. It is a periodic function with respect to v whose period is -. 
Moreover, 



Qi •= hm {v ) — / e^e^ •••e^ (3.1) 

v—>ir}oo ' * 



l<n<i 2 <-<ia<r+l 

is a finite number Thus, we must put 



N_ 

2 
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where the coefficients b^\bj G C are independent of v. Taking note on the 
relation 

m 

lim V — = —5— (3.3) 

n=—m ij 

we can derive the following NLIE for finite Trotter number iV from eqs. ()3.2|) 
and (|2.2()j) by the similar procedures in refs. |Hj and [T3*] . 

f 1 (a) (^) = QS a) 

dy ^ 1} (y + Pi a) - h)Ti a+l \y + Pl a) - 1) 



+ 



+ 



c(a) 2m timr)(v -y- p[ a) )Ti a \y + (3[ a) - 1) 

dy V_Tt% - P[ a) + I )f[ a+1 \y - p[ a) + |) 
c (a) 2vrz tan ^ _ y + pV>)T[ a) {y - p[ a) + i) 

for a G {1, 2, . . . , r}, (3.4) 

where we set fi^ := /^j; the contour (resp. C^) is a counterclockwise 
closed loop around which satisfies the condition y ^ v — (3 < fl l (resp. y 7^ v + 



/3g) and does not surround z{ o) -^+i, -2{j[ a) + ?f, f (resp. z^+^J-i, 

2#' + ^, f) (neZ,fceZ-{0}). 

Next we will take the Trotter limit iV — > 00. We put (v) := limjv-+oo (y) . 
In particular, we have 



T}°\v) = lim f<°>(„) = exp ( S^sinh^) 2 \ 
iv-s-oo \T(cosh 77 — cos(2r]f ))/ 
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Then, we arrive at a system of NLIE, which contains only a finite number of 
unknown functions {T± a \v )}i< a <r : 



i 



(a) 



00 = Qi + 



+ 



£ 
£ 



[ dy yT^jy + f)Tl a+1 \y + f) 
ci a) 2m tan7 j( v -y- *te±H)7} a \y + ^) 
£ dy yT^jy-frf^jy-f) 
2?ri tanr?^ - y + l -^)T} a \y - ^L) 



for a e {1,2,..., r}, (3.6) 



where T^ r+1 \v) = Qx +1 ^ = e T + ; the contour (resp. ) is 



a counterclockwise closed loop around which satisfies the condition y ^ 
v—fl-j^l (resp. y ^ v+f3^) and does not surround z[ flifn+i, — (a+l)i+—, 

f (resp. z[ a) + /?g - i, (o + l)i + f ); /?g = lim^ 0$ = t&i + f; 

z[ a ^ = limjv-^oo^i (« £ ^, H Z - {0})- Note that eq. ()3.6|) reduces to 
Takahashi's NLIE [5 : if r = 1. Although we assumed that q = e v is not a 
root of unity, we think that eq. ()3.6j) will be also valid as it stand even when 
q is a root of unity. One can calculate the free energy per site / by using eq. 
f)3.fij) and the relation 



4 High Temperature Expansion 

In this section, we will calculate the high temperature expansion of the free 
energy eq. ()3.7|) for U g (sl(r + 1)) by our new NLIE eq. ()3.6|) . It is known 
that the high temperature expansion of the free energy for the XXX-model 
was calculated [TT by Takahashi's NLIE up to order 100. 

First of all, we assume the following expansion for large T/\J\: 



where b^(v) = logQi • Substituting eq. (|4.1|) into eq. ()3.6|) . we can obtain 
the coefficients {b„ * 0>)}- We treat r = 1 case separately. 
• r = 1 case: 

In this case, the Hamiltonian eqs. (|2.fi|) and ()2.7j) becomes XXZ chain in 



/ = Jcosh77-Tlog7? i; (0). 



(3.7) 




a e {l,2,...,r}, 



(4.1) 
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magnetic field: 

H + H ch = 2JJ2{SjSJ +1 + S?S] +1 + A(S*S* +1 + -)} 

3=1 

L 

-LQ il +n 2 )-2Q il -vL 2 )Y l S;, (4-2) 

3=1 

where A = cosh?]. The NLIE has especially simple form. Equation ()3.6|) is 

-<■>),„ , .-™ T{°\y-i/2) \ dy 



J c \tanr](v -y -', 



i) tanr7(f — y + i)J T^\y) 2?™ 

(4.3) 

Contour C is a closed loop counterclockwise around 0. To remove trigono- 
metric functions, we use transformation 

v = F(X), y = F(Y) : F(JT) = - taW - zXsinhr ? = ), (4.4) 

By this transformation, points (0, ±z, 00) on v plane move to (0, 00, ±2) on 
X plane, respectively. Then eq. ()4.3j) becomes 

T} X \F(X)) = Q? + QS 2) (1 + X 2 ) I A(X, Y) , 1)( l * Y 

Jc i x (F(Y)) 

Y(1-X 2 Y 2 +2A 2 X 2 (1+Y 2 )) 1 A(l+(l-2A 2 )X 2 y 2 ) 
A(X Y) = V i±Z! + e Vl+A 2 Y 2 

' ^ (i-x 2 r 2 ) 2 -4A 2 x 2 r 2 (i + x 2 )(i + y 2 ) 



j \/i 4- A 2 y 2 

x exp{^(l + F 2 )(A + 6 V1+ y Ay )}, (4.5) 

where we assume that Tl X) {F{Y)) is an even function of Y . New integration 
kernel A is expanded by power series of J/T. Coefficient of Z-th term is 



1 + Y 



2\l 



/! (1 - x 2 r 2 ) 2 - 4A 2 X 2 F 2 (1 + X 2 )(l + Y 2 ) 



Yil -X 2 F 2 + 2A 2 X 2 (1 + F 2 )) //A Vl + ^A 2 ,, /A x/T+yW 
1 + y2 (( + p ) + ( A y 

A(i + (i - 2A 2 )x 2 r 2 ) VTTy^ 

Vi + A 2 F 2 Y ' y Y ' ')' v ' 
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Analytic part of function A at Y = is not relevant in this integral. Then 
we can put 

°° J 1 

A = (analytic part at Y = 0) + -7—7A1, (4.7) 

z — ' l\T l 
1=1 

We get Ai from eq. (J4.6j) as polynomials of X 2 , A 2 and 1/Y; 

A „ 1 + 2A 2 + 2A 2 X 2 

A i = Y> A2 = Y ' 

A 6A + 4A 3 + (3A + 8A 3 )X 2 + 4A 3 X 4 1 + 4A 2 + 2A 2 X 2 

As ~ ys + y ' ~~ y3 

3 + 24A 2 + 8A 4 + (1 + 24A 2 + 24A 4 )X 2 + 8(A 2 + 24A 4 )X 4 + 8A 4 X 6 
+ y ' 

(4.8) 



We can determine b\ 1 \F(X)) successively from the integral equation. These 
are polynomials of X 2 , A and Qf ] /Q? 2 (= 1/(2 cosh £) 2 , h = ( j j 1 -fi 2 )/2 ). 



^(0) = ^-^, (4.9) 
4f)(0)= (l±l^_6A^ i 

q' 1 ' 2 q< 1)4 3q;" 8 

( (j+2A^f)Qf (-1- **-14A«) Qf' 2 

4<0H tf' 2 + 

(40 A 2 + 80 A') Qf' 3 _ 140 A'Qf' 4 

In Table [TJ we give j\of (0) up to j = 10. 

There are a lot of works on high temperature expansions for the XXX 
model; only a few works exist for the XXZ model. Destri and de Vega 's [4*0] 
result contains some misprint. Rojas et al [33] also calculated the coefficients 
up to order 0((^) 4 ). Unfortunately, their formula (eqs. (46)-(48) in ref. 
41J ) contains a misprint of a sign. We note that eqs. ()4.9j) - (j4.11j) agree with 
eq. (46) in ref. [H]. In addition, eq. 1)4.12)1 also agrees with eq. (46) in ref. 



2 A Qf ] 
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Table 1: j\-bf>(0) for j = 1, 10. We put M = /<9i (= 1/(2 cosh £) 2 , 
h = (ti 1 -ti 2 )/2). 

1 2MA, 

2 2Af + 4Af A 2 - 12Af 2 A 2 , 

3 12 M A - 36 Af 2 A + 8 M A 3 - 72 Af 2 A 3 + 160 Af 3 A 3 , 

4 6 M - 36 Af 2 + 48 Af A 2 - 448 Af 2 A 2 + 960 M 3 A 2 + 16 Af A 4 - 336 Af 2 A 4 
+ 1920 M 3 A 4 - 3360 Af 4 A 4 , 

5 60 M A - 960 M 2 A + 2400 Af 3 A + 160 M A 3 - 3600 Af 2 A 3 + 20000 Af 3 A 3 - 33600 M 4 A 3 
+32 Af A 5 - 1440 Af 2 A 5 + 16000 Af 3 A 5 - 67200 Af 4 A 5 + 96768 Af 5 A 5 , 

6 20 M - 600 M 2 + 2400 A'f 3 + 360 Af A 2 - 13320 Af 2 A 2 + 85200 Af 3 A 2 - 151200 Af 4 A 2 
+480 Af A 4 - 23808 Af 2 A 4 + 259200 Af 3 A 4 - 1048320 Af 4 A 4 + 1451520 Af 5 A 4 

+64 Af A 6 - 5952 Af 2 A 6 + 115200 Af 3 A 6 - 873600 Af 4 A 6 + 2903040 Af 5 A 6 
-3548160 Af 6 A 6 , 

7 280 Af A - 18760 Af 2 A + 168000 Af 3 A - 352800 Af 4 A + 1680 Af A 3 
-135408 Af 2 A 3 + 1689856 Af 3 A 3 - 7197120 Af 4 A 3 + 10160640 Af 5 A 3 
+ 1344 Af A 5 - 141120 Af 2 A 5 + 2696960 Af 3 A 5 - 19756800 Af 4 A 5 
+63221760 Af 5 A 5 - 74511360 Af 6 A 5 + 128 Af A 7 - 24192 Af 2 A 7 
+770560 Af 3 A 7 - 9408000 Af 4 A 7 + 54190080 Af 5 A 7 - 149022720 Af 6 A 7 
+ 158146560 Af 7 A 7 , 

8 70 Af - 8820 Af 2 + 117600 Af 3 - 352800 Af 4 + 2240 Af A 2 - 322560 Af 2 A 2 
+5388096 Af 3 A 2 - 26530560 Af 4 A 2 + 40642560 Af 5 A 2 + 6720 Af A 4 
-1145088 Af 2 A 4 + 25079040 Af 3 A 4 - 193277952 Af 4 A 4 + 629959680 Af 5 A 4 
-745113600 Af 6 A 4 + 3584 Af A 6 - 780288 Af 2 A 6 + 24729600 Af 3 A 6 
-292626432 Af 4 A 6 + 1625702400 Af 5 A 6 - 4314562560 Af 6 A 6 
+4428103680 Af 7 A 6 + 256 MA 8 - 97536 Af 2 A 8 + 4945920 Af 3 A 8 
-91445760 Af 4 A 8 + 812851200 Af 5 A 8 - 3775242240 Af 6 A 8 + 8856207360 Af 7 A 8 
-8302694400 Af 8 A 8 , 

9 1260 Af A - 337680 Af 2 A + 8023680 Af 3 A - 50803200 Af 4 A + 91445760 Af 5 A 
+ 13440 Af A 3 - 4072320 Af 2 A 3 + 116998560 Af 3 A 3 - 1030458240 Af 4 A 3 
+3617187840 Af 5 A 3 - 4470681600 Af 6 A 3 + 24192 AfA 5 - 8570880 Af 2 A 5 
+311921280 Af 3 A 5 - 3884186880 Af 4 A 5 + 21973690368 Af 5 A 5 - 58502062080 Af 6 A 5 
+59779399680 Af 7 A 5 + 9216 Af A 7 - 4112640 Af 2 A 7 + 209088000 Af 3 A 7 
-3759436800 Af 4 A 7 + 32286449664 Af 5 A 7 - 144850083840 Af 6 A 7 + 328786698240 Af 7 A 7 
-298896998400 Af 8 A 7 + 512 Af A 9 - 391680 Af 2 A 9 + 30976000 Af 3 A 9 - 835430400 Af 4 A 9 
+ 10762149888 Af 5 A 9 - 75107450880 Af 6 A 9 + 292254842880 Af 7 A 9 

-597793996800 Af 8 A 9 + 501851750400 Af 9 A 9 , 
10 252 Af - 128520 Af 2 + 4445280 Af 3 - 38102400 Af 4 + 91445760 Af 5 + 12600 Af A 2 

-7072800 Af 2 A 2 + 282582720 Af 3 A 2 - 3083330880 Af 4 A 2 + 12395980800 Af 5 A 2 
-16765056000 Af 6 A 2 + 67200 Af A 4 - 42297600 Af 2 A 4 + 1999721280 Af 3 A 4 
-28282381440 Af 4 A 4 + 171577405440 Af 5 A 4 - 475808256000 Af 6 A 4 + 498161664000 Af 7 A 4 
+80640 Af A 6 - 58803840 Af 2 A 6 + 3444222720 Af 3 A 6 - 65241899520 Af 4 A 6 
+570726051840 Af 5 A 6 - 2568513024000 Af 6 A 6 + 5803583385600 Af 7 A 6 
-5230697472000 Af 8 A 6 + 23040 Af A 8 - 20930560 Af 2 A 8 + 1671936000 Af 3 A 8 
-44003635200 Af 4 A 8 + 548674560000 Af 5 A 8 - 3702575923200 Af 6 A 8 
+ 13948526592000 Af 7 A 8 - 27675648000000 Af 8 A 8 + 22583328768000 Af 9 A 8 

+ 1024 Af A 10 - 1569792 Af 2 A 10 + 191078400 Af 3 A 10 - 7333939200 Af 4 A 10 + 131681894400 Af 5 A 10 
-1295901573120 Af 6 A 10 + 7439214182400 Af 7 A 10 - 24908083200000 Af 8 A 10 
+45166657536000 Af 9 A 10 - 34326659727360 Af 10 A 10 
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if we replace the term — 1^/3 3 in eq. (48) with +^P 3 - Moreover, these 
coefficients eqs. (j4.9|) - (j4.12J) reduce to known results (see for example, refs. 
[TT] and 02]) in the rational limit A — > 1. We have plotted specific heat by 
using our high temperature expansion formula of order 99 (see Fig. 
• r > 2 case: 

This case is more involved than r = 1 case since we have to take account 
of the poles in the coefficients {bn\v)}. Thus we need a further non-trivial 
assumption 

&(•)(„) = ( 2 ( sinh ^ 2 V C M(«) (4 13) 

where c$(v) is a polynomial with respect to some trigonometric functions. 
We will present bn\o) for U q (sl(r + 1)) up to order n = 3. 

i?(0) = (4.14) 

^ = (1 + 2^ _ 6A^ + (zi±i|)^, (4 , 5) 
Qr } Qi Q? 



b n m = ( 2A + i f) ^ , (-6 a -12 a 3 ) gf> 2 8 QA 3 gj 2 > 3 
Qf QS 1)4 3Q« 6 

, 8A 3 gj 3) (8A-32A 3 )gj 2) gS 3) (-4A + 8A 3 )Qj 4) 
+ + ms + ^4 • ( 4 - 16 ) 

We note the following relation: bn (0) = ci^(0). We can obtain coefficients 
for r = 2 if we formally set g^ = in eq. f)4.16j) . Note that eqs. ([4.14)1 . 
gUSJ, dHH reduce to eqs. (jO)l . ()4~TU)l . (I4~TT1) if Qj 3) = g^ 4) = 0. Note also 
that eqs. 1)4.14)1 - 1)4.16)1 reduce to eq. (5.6) in ref. ^Hj in the rational limit 
A — > 1. We have plotted the specific heat by using our high temperature 
expansion formula eqs. (J4.14j) - (J4.16j) of order 3 (see Fig. HJ. 



5 Concluding Remarks 

We have derived NLIE with only r unknown functions for the U q (sl(r + 1)) 
Perk-Schultz model. In contrast with traditional TBA equations, our new 
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Figure 3: Temperature dependence of the specific heat C for U q (sl(2)) (fix = 
—fi2 = 0) calculated from Pade approximation (numerator: a polynomial 
of degree 49, denominator: a polynomial of degree 50) for A =2 (thin), 1 
(medium), 0.5 (thick), -0.5 (dashed thick), -1 (dashed medium), -2 (dashed 
thin). We have calculated high-temperature expansion of free energy up to 
99-th order. The case for A = ±1 coincides with ref. [TT] . 



NLIE will be also valid even when q is a root of unity. We also calculated 
the high temperature expansion of the free energy from our new NLIE. 

It will be an important problem to apply our result to physical problem 
such as the thermodynamics of the spin ladder models as in ref. jT3J . In Ref. 
|lb'| . they treated a spin ladder model whose leg part has isotropic interaction 
while rung part has anisotropic interaction. In contrast, our new NLIE has a 
potential applicability to spin ladder models whose leg parts have anisotropic 
interaction. 

We note that our results in this paper reduce to the ones in ref. [T3j in 
the rational limit q = e v — > 1. 



Acknowledgments 

This work is partially supported by Grant-in-Aid for Scientific Research from 
JSPS (no. 16914018). 



15 



1 2 3 4 5 6 

T 



Figure 4: Temperature dependence of the high temperature expansion of the 
specific heat C for U q {sl(A)) (J = 0.48, /ii = 6.75, /i 2 = Hz = Hi = 0): each 
line denotes C for r\ = | (thin),l (dashed thin), (thick), y (dashed thick) 
respectively. The case for r] = was calculated in ref. ^Hj in relation with 
the sw(4) ladder model. 
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